Optical detection of a BCS phase transition in a trapped gas of fermionic atoms 
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o^: 

0\ , Light scattering from a spin-polarized degenerate Fermi gas of trapped ultracold Li atoms is 

,— ( ■ studied. We find that the scattered light contains information which directly reflects the quantum 

-^ ' pair correlation due to the formation of atomic Cooper pairs resulting from a BCS phase transition 

fr\ , to a superfluid state. Evidence for pairing can be observed in both the space and time domains. 

Cn , The realization of Bose-Einstein condensation in trapped atomic gases |l| has generated interest in the atomic 

' physics, quantum optics and condensed matter physics communities. Although the experimental realization of a 

f~| ■ degenerate atomic Fermi gas has not yet been demonstrated, interest in this subject is increasing &-ra]. Of course, 
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O ■ the behavior of a degenerate Fermi gas is remarkably different from a degenerate Bose gas. By analogy vifith the BCS 
Zh theory of superconductivity in metals, it has been predicted that a degenerate Fermi gas can undergo a BCS phase 

transition to an atomic superfluid state if the interatomic interaction in the gas is attractive B . Experiments to trap 
and cool ^Li and ^°K gases into the quantum degenerate regime are underway in several laboratories. 

In this paper, we address the question of how to detect the superfluid state after the BCS phase transition. We 
^ , assume that the Fermi gas has been cooled to near absolute zero, so that all trap levels up to the Fermi energy are 



filled. An attractive interatomic interaction will cause atoms in the vicinity of the Fermi level to form Cooper pairs, 

with each pair composed of two quantum correlated atoms behaving as a new composite Bose particle. These bosons 

automatically undergo Bose-Einstein condensation and form a superfluid. The quantum pair correlation of the Cooper 

'^ , pairs characterizes the superfluid properties of the gas. 

fi ' A promising experimental approach is to prepare a degenerate gas with atoms in an incoherent mixture of two 

internal hyperfine state. Such a mixture allows Cooper pairing via an s-wave interaction, and leads to practically 

attainable BCS-transition temperatures when the scattering length a is large and negative. This occurs naturally for 

^Li H, or can be obtained in the vicinity of a Feshbach resonance for other atoms 0]. We consider here a trapped 

^Li gas in an incoherent mixture of ground states |-|-) = \Ms = 1/2, Mj = 1) and |-) = \Ms = 1/2, Mj = 0) ||]. 

The key to observing the superfluid state is to determine the existence of pair correlations. To achieve this goal, 

{ps ' we propose to use off- resonance light scattering and Fourier imaging techniques. A laser beam with amplitude E^, 

CO , frequency lul, and wave vector k propagating along the z direction is used to illuminate the gas. We take the light to 

CN ' be linearly polarized and tuned near resonance between an S ground state and P excited state. To avoid incoherent 

^i^ , heating of the gas due to spontaneous emission, the magnitude of the laser detuning, S = ujl ~ ojq, is assumed to 

be large. In vector quantum field theory [|9|^l|, the atoms in the light field can be described by a four-component 

atomic field \I/(r) = -0+1+) +-0-1—) +^e+|e-l-)-l-'0e-|e— ) with ip± denoting atoms in the ground-state hyperfine levels 

C^ ' |±), and ipe± in the corresponding excited-state hyperfine levels. For large S, the excitcd-statc components can be 

G , adiabatically eliminated, yielding a total atomic polarization operator with positive-frequency part B 



^ '■ P^+) (r, t) = -p^:^p{r, i)e— *, (1) 

O- , , 

O ^ where /5(r,t) = tp'^{r,t)'ip^{r,t) + ^_{r,t)^_{r,t) denotes the total atomic density operator in the ground state, p 
^ I the matrix element of the atomic dipole moment, and r a location in the gas. Light propagation is determined by the 
: '"j I atomic polarization operator (nl) and the wave equation 

..=^: V^eW-?^^=Mo^^. (2) 

The solution to Eq. (g) can be expressed as 

E(+)(R,t) ==El;+'(R,i)e-*"^* + E^+^e*'''^-''^^*, (3) 

where Eg (R, t) is the scattered field at position R. For R = |R| 3> |r|, the scattered field has the form [p|,pT[ 

E^+) (R, f ) = fc2 £^ / d3^e-''=^-^ [P(+) (r, i) - R • p(+) (r, t)Rj , (4) 



where the directional unit vector R = H/R. From Eqs. (|]) and (^, we see that the scattered field depends on the 
density operator of the gas, so that the averaged spectral intensity of the scattered field received by a photodetector 
contains the second-order correlation of the atomic field operators ||l^ 

(p(r, t)p(r', t')) « (p(r, t)){p{r\ t')) + G(r, r', t, i'), (5) 

where "(. . .)" denotes the quantum mechanical expectation value. The first term in Eq. dq), which depends on the 
total averaged density, describes the contribution to the scattered field by the normal ground-state component. The 
second term, 

G(r,r',i,t') = -2(^_(r,t)7A+(r',t'))(V'l(r,t)V'Ur',t')), (6) 

gives the quantum pair correlation function arising from the formation of Cooper pairs in the superfluid state. 

The contribution of the laser field E^ in Eq. (Q) can be removed by imaging the cloud with a dark ground technique, 
as discussed in Refs. U%- If a plane located a distance zq from the atoms is observed in this way, the spectral and 



spatial intensity distribution measured on the detector will be |13 



/oo -1 pT 

dre^'^'^y di(E^"HRo,i)-E^+)(Ro,t-FT)), (7) 

where 2T is the time interval used for detection, and Rq = (R^, zq) is a point in the image plane. Equation (H), along 
with relations (ll|) and dsl) , gives the spatial-temporal correlation function of the light field 

(E^-^(Ro,t)-E^+^(Ro,t + r)) = J^^';^^)^ [/i(R±, t, r) + /2(R^, ^,t)] e''" ^-, (8) 

where II = E^ ' ■ E^ ' is the intensity of the incident light and 7 is the natural linewidth of the transition. The 
functions /i and I2 are defined as 

/i(Rj.,i,T)=y'y'd2j.^dV^e-''=^-(-^--^^)/^«(p(r^,i))(p(rY,t + r)), (9) 

and 

l2{-R±,t,T)^ J d^^e-^'^'^^^/'" J d^r^G{r^,r^-^,t,t + T), (10) 

where the relative distance between atoms is denoted by ^ = r^ — r'j^ . The function /i describes the signal from the 
normal component of the gas and I2 the signal from the Cooper pairs. In general, I2 is much weaker than /i since 
the averaged density of atoms in the normal component is far larger than that of the pairs. 

The averaged density and the quantum pair correlation function can be found using vector quantum field theory 
P|. In the off-resonant light field, the degenerate Fermi gas is described by the coupled quantum field equations 






ih-^ = [Ho -ii+ + Vl- inT/2)^j+ - A(r)^^ 

-{Ho- ^i. + VL + inT 12)^1 -^{v)^+, (11) 



where Hq = — ^^V^ -I- \muj'^r'^ is the free Hamiltonian of the trapped Fermi gas, Vl = ?iri^/4(5 is the light-induced 
potential, F = "^VP /AS"^ is the rate for spontaneous emission, ^± are the chemical potentials of the two internal states, 
and A(r) = (47r|a|?i^/m) (■)/'- (r)'0+(r)) is the BCS energy gap function |^. The Rabi frequency of the light field is 
f2 = |p • Ei/?i|. For simplicity we consider the simple case /i+ = /i_ for equal number of atoms in each spin state and 
introduce the renormalized chemical potential ^ = /i+ — Vl- Further, a large laser detuning and a weak intensity allow 
F <C /i, A so that destruction of Cooper pairs by spontaneous emission and interactions involving excited-state atoms 
PI can be neglected. Employing an approach similar to that adopted in BCS theory, we approximate the solutions of 
Eqs. (|ri|) by 

V'±(r, t)^Y. («n(r)S„±e-^^"*/'' ± v^{v)hl^e^''""^) , (12) 



where &n± are generalized Bogoliubov quasi-particle operators and En the excitation energy for the mode indexed 
by n. The superfluid state of the degenerate Fermi gas is characterized by the BCS ground state |$_bcs) with the 
property 6n±!^BCs) — 0. From Eqs. ( |ll| ) with the dissipative terms ignored, the transformation coefficients {un,Wn} 
satisfy the celebrated Bogoliubov equations 

{Ho - /i)un(r) + A(r)i;„(r) = EnUn{v) 
-{Ho - /i)vn(r) + A(r)w„(r) = EnVn{r). (13) 

The total averaged density can be expressed as {p{r,t)) = {^bcs\'4'+'>JJ+ + V'-V'-I'i'-Bcs) — 2X]n l'^n(r)P and the 
quantum pair function is 

G{r,r',t,t') = 2^w„(r)t;„(r')^m(r')z;m(r)e-*(^"+^'")(*-*')/''. (14) 

nm 

The average density and pair function can be calculated by self-consistently solving Eqs. (|l^). In the normal 
degenerate ground state, energy levels below the Fermi level Ep are occupied, while those above are empty. The 
effect of interatomic interactions is to cause scattering between nearby energy levels, which creates an energy shell 
near Ep where normally unoccupied states in the normal ground state acquire an amplitude to be occupied, and 
states below Ep have some amplitude to be unoccupied. The stronger the interatomic interaction is, the wider the 
energy shell and the more atoms are available to form Cooper pairs. Physically, the coefficients Un and v^ in Eqs. ( |l^ ) 
determine the amplitudes for atoms to be scattered into the pair states. To evaluate these amplitudes, we expand the 
coefficients as u^ = X]q Unq't'q and Vn = J2q "nq^q, in terms of the eigenstates (/>q of the single-atom Hamiltonian Hq. 
In principle the sum over q in the coefficients should extend from zero to infinity. However, we should note that in 
the BCS theory |§,^, Equation (13) is a direct result of the Born approximation by replacing the realistic non-local 
interatomic interaction V{r) by a local contact potential V{f) = inh a6{r)/m. We know that the Born approximation 
is only valid for low-energy scattering. The invalidity of the approximation in high-energy scattering regime produces 
an ultra-violet divergence in the BCS theory. In the case of superconductivity, the ultra-violet divergence naturally 
vanishes by considering the fact that the phonon-exchange induced interaction between electrons can be cut-off in 
the Debye frequency. However, in the case of degenerate Fermi gas of atoms, to avoid the ultra-violet divergence, an 
exact theory for superfluid phase transition must take the realistic shape of the exact non-local triplet potential into 
account. Recently two independent approaches to remove the ultra-violet divergence in the BCS theory of degenerate 
Fermi gas of atoms have been proposed P,0| . One is to renormalize the interaction potential in term of the Lippman- 
Schwinger equation |2] and the other is to employ the more exact pseudo-potential approximation 0]. However for 
a first guess, the Born approximation provides a simple and reasonable way to evaluate the gap energy and the pair 
correlation if an appropriate momentum cut-off is introduced to remove the ultra-violet divergence. Now the question 
is how to choose a physically valid momentum cut-off fikc- To determine the cut-off range, we must use the fact that 
Born approximation only gives the correct evaluation in the low-energy scattering regime with k\a\ < 1. Hence the 
validity of the present theory based on Born approximation requires a cut-off fee < 1^1 • For ^Li atom, this is in 
the order of Fermi wave number kp. With such a cut-off, we numerically evaluate the energy gap, the total averaged 
atomic density and the quantum pair function. 

To be concrete, we assume that N = 2 x 10^ ®Li atoms in each spin state are confined by a magnetic trap with an 
oscillation frequency a; := 27r x 150 Hz. With these values, Ep « lOOhu; « 740 nK, and the peak value of the energy 
gap is A(0) « 5tiuj = 36 nK. For a degenerate Fermi gas in a harmonic trap, the characteristic size of the average 
density is given by the Fermi radius rp = [2Ep/muj'^]^/^ « 48 fira |1J], while the length scale of the pair correlation 
function is Tc ~ kp^ , where kp = {2mEp/h )^'^ « 27r x 6800 cm^^ is the Fermi wavenumber. The numerical result 
for the correlation function is shown in Fig. n, along with the spatial variation of the energy gap. 

We need emphasize that in the homogeneous gas, the correlation length (pair size) at zero temperature is defined 
in terms of the so-called coherence length ^c = hvp/A{0). The coherence length determines the region where the pair 
function extends M. However within the region, the pair function still contains shorter oscillation structure which 
has the scale re. In fact from our numerical result for the trapped gas, we see that the pair function indeed varys 
with such a length scale. Now we will explain how such a scale can be observed by optical imaging. 

Assuming that a plane at zq = 2 cm is imaged with unit magnification and with a transition wavelength A = 670 
nm, the image size is zo/krp '^ 0.09 mm for the normal component and zo/kvc '^ 1.9 cm for the pair component, 
differing by a factor of 2Ep/huj. The calculated images for a gas below and above the critical temperature for the BCS 
phase transition are shown in Fig. 0(a) and (b), respectively. It is seen that when the transition occurs, a spatially 
broadened image appears. The physical situation is depicted in Fig. 0, where the small-scale structure induced by 
pairing causes light to scatter at a larger angle than that scattering from the cloud itself. 



The normal signal is produced by coherent scattering and is therefore proportional to {2N)'^, as can be verified 
by reference to Eq. (9). The pair signal, however, arises from spontaneous Raman scattering between pairs above 
and below the energy gap, and is found using Eq. (10) to be proportional to the number of pairs Np. This number 
is determined by the number of atoms in an energy shell of width A centered on Ep, so Np w 3NA/Ep. For the 
parameters given above, Np « 3 x 10'' and the ratio of the peak signal intensities is /2(0)//i(0) w 2 x 10~^. It is 
difficult to experimentally measure a signal with such a large dynamic range, but the pair signal can be revealed 
by using a nearly opaque spatial filter to attenuate the normal signal. If the diameter of the filter is chosen to be 
approximately equal to the spatial dimension of the normal signal image, it will affect only the central region of the 
pair signal, and both contributions can be observed with the same intensity scale. 

Finally, we calculate the scattered light spectrum. For the normal degenerate ground state, a single spectral line 
is obtained at the frequency of the incident light. For the superfluid state, the spectrum exhibits a double-peaked 
structure as shown in Fig. 0. The coherent peak is from scattering by the normal component. The frequency shift 
of the sideband line is approximately twice the gap energy, confirming that the sideband is due to Raman scattering 
by pairs. The long oscillating tail of the sideband is due to modulated broadening from the center of mass motion of 
atoms at the trap frequency. Hence, the presence of the shifted peak provides another effective method to detect the 
BCS phase transition and can be used to directly determine the gap energy. 

The theory presented here was simplified by the neglect of spontaneous emission, permitting, for example, the 
assumption that A remains constant during probing. However, the pair signal depends on breaking pairs by incoherent 
spontaneous Raman scattering, and thus requires spontaneous emission. The theory is therefore valid only in the 
weak-signal limit, where F ^ T~^. Larger signals could be obtained experimentally by allowing F ~ T~^, but 
quantitative interpretation would then be more difficult. 

In conclusion, we have studied off-resonance light scattering by a trapped degenerate Fermi gas. The results show 
that both spatial imaging and the scattered light spectrum give clear signatures for the BCS phase transition to a 
gaseous superfluid state. 
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FIG. 1. The normalized spatial distributions of the equal-time quantum pair correlation function, with the distance ^ 
between atoms scaled by the Fermi wavenumber kp- The dashed curve is the correlation due to Cooper pairs with center of 
mass located at the trap center and the solid curve is the average contribution of all Cooper pairs. The inset shows the spatial 
dependence of the energy gap, normalized to the center of the trap and with position scaled by the Fermi radius rp- 





FIG. 2. The spatial image measured a distance zq from the atoms, with (a) the trapped Fermi gas in the normal degenerate 
ground state and (b) in the superfluid state after the BCS phase transition. In both images, the central peak is clipped and 
actually extends by a factor of ~10^ above the axes shown. The radial size of the normal component is approximately zo/krp, 
while the pair component is larger, extending to zo/krc. 





FIG. 3. Schematic of the imaging technique. The white area represents the incident probe laser, the dark gray the hght 
coherently scattered by the cloud, and the light gray the light scattered by Cooper pairs. The small length scale of the 
pair structure scatters light at a relatively large angle, so by measuring the intensity in the far field, the components can be 
distinguished. Dark ground imaging techniques can be used to eliminate the contribution of the probe laser itself [[L2| . 
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FIG. 4. The normalized scattered light spectrum of the scattered field. The frequency shift Av is scaled by the trap 
frequency co. For the solid curve, the gap energy A ~ 5 hu; — 36 nK, and for the dotted curve A = 72 nK. A spatial filter with 



a transmission of ' 
component. 



TO is used to reduce the strong signal from the normal component to the same level as that from the pair 



